Shape Deformation driven Structural Transitions in Quantum Hall Skyrmions 
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The Quantum Hall ground state away from v — 1 can be 
described by a collection of interacting skyrmions. We show 
within the context of a nonlinear sigma model, that the clas- 
sical ground state away from v — 1 is a skyrmion crystal with 
a generalized Neel order. We show that as a function of fill- 
ing v, the skyrmion crystal undergoes a triangle — > square 
— » triangle transition at zero temperature. We argue that 
this structural transition, driven by a change in the shape of 
the individual skyrmions, is stable to thermal and quantum 
fluctuations and may be probed experimentally. 



PACS: 73.40.HK, 67.57.Fg, 64.70.Kb, 75.10.Hk, 67.80.-s 



Following a recent suggestion J 1 - 3] , there has been 
growing experimental evidence p- 3 [that the charged 
quasiparticle excitations about the v = 1 quantum 
Hall state in GaAs heterostructures are extended objects 
called skyrmions with spin significantly greater than 1/2. 
Recent OPNMR Q and optical magneto-absorption g] 
experiments observe a sharp fall in the spin polarization 
as the filling factor is changed from v = 1, consistent 
with the existence of skyrmions. 

What is the ground state of a system of interacting 
skyrmions in two dimensions ? Specific heat measure- 
ments JTJ] on GaAs heterojunctions at v = 0.77, show a 
sharp peak at a temperature T c rs 30mK. It has been 
suggested M that this anomaly may be associated with 
the freezing of skyrmions into a crystal lattice. Previous 
theoretical work JqM have analysed the crystalline state 
of skyrmions. Using a mean-field analysis of electrons 
confined to the lowest Landau level, Fertig et. al. M 
claim that at T = 0, the skyrmions form a square lattice 
with a Neel orientation ordering. On the other hand, by a 
mapping onto an effective nonlinear sigma model, Green 
et. al. |9[] conclude that the skyrmions form a triangular 
lattice with a generalized Neel ordering. 

In this paper, we calculate the T = classical phase 
diagram of a system of interacting skyrmions using an 
effective classical 0(3) nonlinear sigma model (NLSM). 
Within a variational scheme, we show that at values of 
v sufficiently away from 1, the ground state is a trian- 
gular lattice of skyrmions with a three sublattice gener- 
alized Neel orientation order. As we approach v = 1, 
we find a structural transition to the square lattice with 
Neel orientation order. This structural transition is ac- 
companied by a change in the shape of the individual 
skyrmions and a jump in the spin polarization. In the 



dilute limit, v — > 1, we find a reentrant triangular phase 
which is different from the previously encountered trian- 
gular phase! We also provide arguments for the stability 
of these phases to thermal and quantum fluctuations. 

The low energy excitations about the ferromagnetic 
v = 1 ground state in GaAs hetrostructures are described 
by an NLSM plfl , in terms of a unit vector field that rep- 
resents the local spin polarization. Neutral excitations 
correspond to spin waves, while charged excitations cor- 
respond to defect configurations (skyrmions) with topo- 
logical charge Q (the electronic charge density is equal to 
the topological charge density 47rg(x) = n • (d x n x d y n) 
in the long wavelength limit [Q ) . The topological charge 
Q = J (j'(x) d 2 x is always an integer and counts the num- 
ber of times the spin configuration n(x) wraps around 
the unit sphere. 

For convenience we work in a notation where the 
unit vector field n(x, y) is replaced by a complex field 
w(z = x + iy,z = X — iy), obtained by the stereographic 
projection of the unit sphere onto the complex plane. 
Thus w = cot(|) e 1 ^, where 8 and <j> are the polar angles 
of the unit vector n. For static, classical spin configura- 
tions, the energy functional may be written as, 

E [w, w] = E [w, w] + E z [w, w] + E cou i [w, w] . (1) 



The first term on the right 

E [w,w] = 7 



-(dzwd z w + d z wdzw) (2) 



(1 + ww)' 



describes spin exchange. The Zeeman coupling of the 
spins to the external magnetic field B is given by 



E z [w, w] 
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(3) 



while the charged quasiparticles interact via a long-range 
coulomb interaction, 



E C oul [w, w] 



q(z,z)q(z',z') 



(4) 
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The energy functional Eq. (fil) has been expressed in di- 
mensionless variables, with all energies in units of the 
cyclotron energy hio c , uj c = eB/m*c (m* is the elec- 
tron band mass) and all lengths in units of the magnetic 
length l c , 2irl 2 — hc/eB. The parameters in the energy 
functional are 7 = e*/(16v27r), g* = gusB /fk>J c and 
e* = (e 2 /Kl c )(l/hu! c ) (K is the dielectric constant of 
GaAs). 



In the absence of Coulomb and Zeeman interactions, 
any (anti)meromorphic function w(z) is a solution of the 
resulting Euler Lagrange equations [jll| . The topological 
charge is simply Q — J^ . n% , where i runs over the poles of 
w and n% is the degree of the i th pole. The one skyrmion 
solution given by 



w(z, z) 



\e l 



£ 



(5) 



clearly has a Q = 1. The spin and charge distributions 
are centred at £ and fall off as power laws with a scale set 
by A. The XY component of the spin at x is oriented at 
an angle f2 to the position vector x. The Z component 
of the spin Sz, however, diverges logarithmically. 

In the presence of the Zeeman and Coulomb interac- 
tions, Eq. (0) is no longer the minimum energy solu- 
tion. These terms destroy scale invariance, and generate 
a 'size' for the optimal skyrmion, leading to a finite Sz- 
It is natural to try the following variational form for the 
single skyrmion solution which minimises Eq. (hi), 



w(z, z) 
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(6) 



With A and k as variational parameters, we have checked 
that the lowest energy single skyrmion solution has a 
k > 0. The topological charge remains 1, since the con- 
figuration in Eq. (0) can be smoothly deformed to the 
configuration in Eq. (|j|). 

A system of N identical Q — 1 skyrmions centred at 
{£/} with orientations {£!/}, can now be parametrized by 



N 



w(z,z) = Y^ 



Xe 
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(7) 



7=0 



The topogical charge Q = N, however we observe that 
even for k = 0, the spin Sz is finite. This illustrates 
the fact that the Z component of the total spin of over- 
lapping skyrmions is not the sum of the individual spins. 
Our crystalline ansatz corresponds to placing the {£/} on 
a triangular or a square lattice. We have studied both 
the ferro-oriented (SI/ = 0, V7) and generalized Neel ori- 
ented (described below) configurations. We find that the 
ferro-oriented configurations always have a higher energy 
and so shall ignore them for the rest of this paper. Since 
the square lattice is bipartite, the Neel configuration is 
characterised by $1/ = for the A sublattice and 0/ = 7r 
for the B sublattice. Likewise, for the tripartite triangu- 
lar lattice, the generalized Neel ordering is obtained by 
assigning SI/ = 0, 27r/3, 47r/3 to the A, B and C sublat- 
tices respectively. These orientation assignments emerge 
naturally in the k = limit of Eq. (m) , when it reduces to 
an elliptic function. The sum of the residues of an elliptic 
function in the unit cell should equal zero. In the present 
context this is precisely the generalized Neel condition 
that the sum of the orientations of all the skyrmions in 
a unit cell is zero (mod 27r). 



To find the classical ground states, we numerically 
compute the energy (accurate to 1 part in 10 6 ) using 
our crystalline ansatz and minimize with respect to the 
variational parameters k and A for a given lattice and 
orientational ordering. We have chosen typical values 
of the carrier concentration (n — 1.5 x 10 11 cm -2 ) and 
magnetic fields (B = 6 — 10 Tesla) and have varied v 
by tilting the magnetic field with respect to the normal 
keeping its magnitude fixed. The lattice spacing a for a 
given lattice type is fixed by the value of v (e.g., for a 
square lattice, a = a/27t/|1 — u\) . 
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Fig. 1 The minimised energy E(k) as a function of k 
for the triangular lattice at v — 1.15, v = 1.04 and v = 1.01. 
The magnetic field B=TT. 

To highlight the role of the 'size' parameter k, we refer 
to Fig. 1, a plot of E(k), the energy of the Neel triangu- 
lar lattice minimized with respect to A, as a function of 
K for three different values of v. As v — ► 1, the lattice 
gets very dilute, and the energy E(k) attains its mini- 
mum at a value of k close to the single skyrmion value. 
Next consider the dense region when k _1 > a (strongly 
overlapping skyrmions). The nonadditivity of spins in 
this region, implies that Ez is insensitive to k. The self- 
part of E cou i and Eq however, increase linearly with in- 
creasing k (for small k) as in Fig. 1. Keeping a fixed, 
a further increase in k will separate the skyrmions from 
each other, giving rise to another minimum at the single 
skyrmion value of k. These minima become degenerate 
at v k, 1.05. A further decrease in v makes the energy of 
the second minimum lower. This suggests that there are 
two different triangular crystals associated with k = 
("fat" skyrmions) and K > ("thin" skyrmions). For the 
square lattice, the "thin" skyrmion solutions have lower 
energy throughout the v range of physical interest. 

Figure 2 shows the minimised energies of the Neel 
square and the Neel triangular crystals as a function of 
v. When v is away from 1, a triangular lattice with n = 
(Ai phase) has the lowest energy. The smaller Zeeman 
energy of the competing square lattice with k > is off- 



set by Eq and E cou i. As v approaches 1, E cou i decreases, 
resulting in a weak first-order transition (slope disconti- 
nuity) to a square lattice with k > at \v — 1| ss .05 (for 
aB = 7T). Thus the structural transition between the 
Ai and the square phase is accompanied by a change in 
the size of the individual skyrmions. In the region of pa- 
rameter space which favours the square phase, the k > 
triangular lattice loses out on a higher Eq and Ez- As 
\v — 1 1 gets very close to zero, the difference in energy be- 
tween these two lattices vanishes. In the limit of extreme 
dilution, v — » 1, the energy of the A-skyrmion configu- 
ration can be evaluated as an expansion in (ko) , 



E 



N 



NEx 
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(8) 



where E\ is the single skyrmion energy and the coulomb 
interaction is between point charges. To leading order, 
E^ is clearly minimised by placing the charges on a tri- 
angular lattice (A2 phase), however it could have any ori- 
entational order. Higher order terms favour a Neel orien- 
tation. Thus a second structural transformation should 
occur between the square Neel lattice and a triangular 
Neel lattice with k > (A2 phase) via a weak first order 
phase transition as v approaches 1. 
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Fig. 2 Energy density (in units of meV/ l c ) as a function 
of v for the triangular (Ai) and square crystals at B — 7 T. 

In the vicinity of the Ai-square transition, we find that 
k _1 ~ a. A simple scaling argument shows that the 
shape stiffness of a single skyrmion is of the same order 
as the elastic stiffness of the crystal in the neighbourhood 
of the transition. The single skyrmion energy is of the 
form 



Ei = 7C1 



g*c 2 R 2 



e c 3 



(9) 



where R sets the scale for the size of the skyrmion 
and ci,C2 and a are constants of Oil). At R ~ a, 
the shape stiffness, d 2 Ei/dR 2 \ R=R , (where R minimises 



i?i) is comparable to the elastic stiffness of the crystal 
~ e 2 /Ka 3 . The novel feature of this structural transi- 
tion is that it is caused by the shape deformability of the 
"atoms" , revealing a richer physics than that of Wigncr 
crystallisation. 

A computation of the spin polarization of the skyrmion 
crystals at T = (Fig. 3), agrees well with the experi- 
mental data of p| obtained at I.5K. Since the spin po- 
larization is a measure of the size R of the individual 
skyrmions, this indicates that R does not change signif- 
icantly over this temperature range. An interesting fea- 
ture is that the Ai (k = 0) to square (k > 0) structural 
transition is accompanied by a discontinuity of about 5- 
10% in the spin polarization, and so may be probed by 
accurate spin polarization measurements. Since thermal 
fluctuations and the presence of quenched disorder smear 
out this discontinuity, it may be necessary to go to very 
low temperatures to see this effect. 
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Fig. 3 Spin polarization as a function of v at B — 7 T 
for the square and triangular lattices compared to the exper- 
imental points taken from Ref. W. The structural transition 
predicted by our theory is accompanied by a discontinuity in 
the spin polarization at v ~ 0.95, as seen from the dashed 
curve. 

How 'good' is our classical variational ansatz Eq. R ? A 
more general ansatz for the individual skyrmions would 
involve deviations of the shape from circular symmetry, 
and giving each skyrmion an arbitrary topological charge 
Q. It is easy to see that both these generalisations lead 
to an increased energy, in the region where the skyrmions 
do not overlap. 

How do quantum and thermal fluctuations affect the 
classical T — phase diagram presented above? The 
qualitative features of the phase diagram in the T — v 
plane (Fig. 4) may be glimpsed from general arguments. 
Moving away from v = 1 along the T = axis, re- 
duces the lattice spacing, and shrinks the skyrmion size 
(due to coulomb repulsion). At a critical lattice spacing, 
zero-point fluctuations would melt the crystal. Quantum 



melting in a strong magnetic field would ensue when the 
magnetic length l c ~ aa, where a ~ 0.1 — 0.2. This 
corresponds to an energy scale of (7i/0.1a) 2 m _1 , where 
m is the mass of the skyrmion. At the melting transi- 
tion, this should compare with the coulomb energy e 2 /a. 
We estimate the quantum melting transition (QM) to 
occur at v ~ 0.8. This melting into a quantum oriented 
liquid (QL) is most likely continuous with an intermedi- 
ate quantum hexatic phase. Beyond this ofcourse, the 
description in terms of skyrmions breaks down as v ap- 
proaches the next quantum hall plateau, e.g. 2/3. In 
addition, we encounter a new quantum orientation dis- 
order transition (QOD) in the limit of extreme dilution. 
The oriented (Neel) crystal is characterised by power-law 
correlations in the sublattice orientation. Quantum fluc- 
tuations would destroy this order when the fluctuations 
in the orientation become of the order of 1-k. i.e. when 
h J II ~ J (a), where I is the moment of inertia of the 
skyrmion and J (a) is the energy cost in changing the ori- 
entation (and is a decreasing function of a). This leads 
to a quantum disoriented crystal in the dilute limit via a 
Kosterlitz - Thouless transition. 

An increase in temperature, T, would result in an ori- 
ented crystal with quasi-long range order in position and 
orientation. The advantage that the square phase had 
over the A2 phase at T — now diminishes, since the 
renormalized 7 gets weaker and the coulomb interaction 
remains relatively unaffected. Thus at higher tempera- 
tures, the square phase disappears, giving rise to an iso- 
structural (first-order) phase boundary, where the two 
triangular Neel phases, Ai and A2 meet. This first-order 
line terminates on a continuous melting curve (a Wigner 
crystal estimate || gives T m ~ a/|1 — f |) which occurs 
via a defect mediated mechanism (KTNHY [Q ) . In this 
case a hexatic phase (H) intervenes between the solid 
and the classical liquid (CL) phase. Indeed there has 
been a recent suggestion [|13J that the presence of an iso- 
structural critical point might reveal a hexatic phase in 
the vicinity. 
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Fig. 4 Conjectured Phase diagram as a function of filling 
v and temperature. The solid and the dashed lines repre- 
sent first-order and continuous boundaries respectively. The 
various phases shown are defined in the text. 



The application of hydrostatic pressure constitutes an 
additional 'field variable' H]. Pressure induces wave- 
function overlap which leads to a reduction in g, thus 
reducing the Zeeman contribution. We find that when 
g — » 0, the square Neel phase ceases to be a ground state. 

As in the case of the Wigner crystal at v ~ 0.2 p5[ , 
we expect that disorder in the GaAs will not destroy 
the (quasi)-long range crystalline order of the skyrmions. 
However, the weak first-order structural transitions re- 
ported above, will be rendered continous. We note that 
the length scales a and kT 1 are both around 10 3 times 
larger than the GaAs lattice spacing, making the sub- 
strate essentially a continuum. Sliding of the skyrmion 
crystal would however be prevented by pinning to the 
ever present disorder in GaAs. 

We would like to thank B. I. Halperin for discussions. 
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